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Abstract We study small clusters of bosons, A = 2,3, 4, 5, 6, characterized by a reso- 
nant interaction. Firstly, we use a soft-gaussian interaction that reproduces the values of 
the dimer binding energy and the atom-atom scattering length obtained with LM2M2 
potential, a widely used 4 He- 4 He interaction. We change the intensity of the potential 
to explore the clusters' spectra in different regions with large positive and large neg- 
ative values of the two-body scattering length and we report the clusters' energies on 
Efimov plot, which makes the scale invariance explicit. Secondly, we repeat our calcu- 
lation adding a repulsive three-body force to reproduce the trimer binding energy. In 
all the region explored, we have found that these systems present two states, one deep 
and one shallow close to the A—l threshold, and scale invariance has been investigated 
for these states. The calculations are performed by means of Hyperspherical Harmonics 
basis set. 
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1 Introduction 



Systems composed by few particles having large value of the two-body scattering length, 
a, with respect to the natural length, I, fixed by the inter-particle potential, have 
been the object of an intense investigation both from a theoretical and experimental 
point of view (for recent reviews see Refs. [TJ[2 3 ). The interest is driven by their 
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universal properties; the behavior of observables do not depend on the microscopical 
characteristics, namely the inter-particle potential, but only on symmetries. 

In the limit a/l —> oo, known either as resonant a — > oo or as scaling £ — > limit, 
the few-particle systems display discrete-scale symmetry as brought out by Efimov in 
his works in the early 70's [HE]. The scaling factor is usually written as A = e^^ s ° , and 
for three-identical bosons sq = 1.00624 and A « 22.69. This symmetry implies that 
all observables can be written as an universal log-period function of the dimensionless 
variable <m*, where k* is a three-body parameter encoding the high-energy (short 
distance) physics, which enters, at leading order, only through this parameter. The 
fact that the observables' behavior is governed by discrete-scale symmetry is know 
as Efimov physics. For ^4 = 3, these properties have been studied for large positive 
and large negative values of the scattering length in the (a _1 ,fc) plane, with k = 
sign(_E)[|_E|/(fi 2 /m)] 1 / 2 , constructing the Efimov plot [5]. This plot is useful in the 
identification of discrete-scale symmetry; in fact, if one introduce the radial, H, and 
angular, £ Efimov variables by k = i/sin^ and a -1 = ffcos£, the scale invariance 
reads H — > \~ 1 H, and £ — > £. 

In the present work we extend our previous analysis of the A = 4 — 6 bosonic 
spectrum [6] to the (a _1 ,fc) plane. We have modified the strength of the LM2M2 [7] 
potential in order to cover the region of negative values of a up to a°_, with this value 
indicating the threshold of having a three-body system bound. We have also increased 
the intensity of the interaction in order to extend the analysis to positive values of a 
in which the universal character of the system starts to be questionable, i.e, when the 
ground-state E® approaches the natural energy Eg = —h 2 /ml 2 , which delimits the 
Efimov window. 

We used the LM2M2 potential to fix the two-body soft-core potential as in discussed 
in Refs. [HUB]; this has been possible because of the scale separation between the 4 He- 
He scattering length, a = 189.41 a.u., and the natural length £ = 10.2 a.u., which is the 
van der Waals length calculated for the LM2M2 potential [I] . In the three-body sector, a 
three-body soft-core potential is required to reproduce the ground-state-binding energy 
of the helium trimer given by the LM2M2 potential. 

We have performed our numerical calculations in systems with A > 3, by means 
of the non-symmetrized hyperspherical harmonic (NSHH) expansion method with the 
technique recently developed by the authors in Refs. 9 10 11 12 . In this approach, 
the authors have used the Hyperspherical Harmonic (HH) basis, without a previous 
symmetrization procedure, and on the representation of the Hamiltonian matrix, as a 
sum of products of sparse matrices, well suited for a numerical implementation. 

As a result, we have observed that in all the region explored the A — 4, 5, 6 systems 
present two states, one deep and one shallow close to the Ej L _ 1 threshold. This analysis 
confirms, at least in one zone of the Efimov plot, previous observations that each Efimov 
state in the ^4 = 3 system produces two bound states in the A = 4 system, and extends 
this observation to the A = 5, 6 systems. 

2 Potentials 

In our calculation we used h 2 /m = 43.281307 (a.u.) 2 K as mass parameter. The LM2M2 
interaction has been modified in the following way 



V x (r) = A • Vlm2M2(0 , 



(1) 
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and we have varied A from A = 0.883, where a = a_ = —43.84 a.u., up to A = 1.1 
corresponding to a = 44.79 a.u. The unitary limit is produced for A w 0.9743. 
The two-body gaussian (TBG) potential is 

V(r) = V e~ r2/R o , (2) 

with range Rq — 10 a.u., and we have varied the strength Vq in order to reproduce the 
values of a given by V\(r). In the three-body sector we need an hypercentral-three- 
body (H3B) interaction to better describe the A — 3 system obtained with the modified 
LM2M2 potential 

W( P1 23) = W e -"?23/Po ( ( 3 ) 

with the strength Wq tuned to reproduce the trimer energy E® obtained using V\{r). 
Here p| 2 3 — § ( r i2 + r 23 + r 3l) ^ s the hyperradius of three particles and the range of 
the three-body force po = ^0 

3 Results 

We have solved the A = 3 problem for bound states using the modified LM2M2 po- 
tential, and then we used the resulting energies to fix the strength of the H3B force. 
Than we have diagonalized the Hamiltonian for A = 3, 4, 5, 6 bodies using the TBG 
and TBG+H3B potentials. The results are given in Fig. [T]in two (a~ 1 ,k) plots, which 
have been scaled to shrink the scale factor to y/X « 4.8. 

In the case with only the TBG potential, upper panel of Fig. \T\ we observe that 
the spectrum of the systems A = 4, 5, 6 presents two bound states, one deep and one 
shallow, for all values of a studied. These calculations confirms the prediction for A = 4 
of a pair of tetramers attached to a trimer [1311141115] . and extends the observation to 
A = 5 and A — 6 systems. 

When the repulsive three-body force is included, lower panel of Fig. [TJ the spectrum 
moves up and we can observe that the excited state E\ disappears for A = 5, 6 for 
negative values of the scattering length as a approaches a_. This behaviour is sensitive 
to the range of the three-body force po, and it has been deeply investigated in Ref. [16] . 

In Fig. [2] we investigate the scale invariance using Efimov-polar coordinates: for a 
fixed value of £ we report the ratio (E\/ E^) 1 ^ 2 calculated with TBG potential. For 
negative values of a, the ratios tend to be constants in agreement with discrete-scale 
invariance; we note that, even if the ratio for A — 3 tends to be constant, the value is 
lower than the expected A. For a > 0, which corresponds to the shaded zone in Fig. [2] 
the non-universal behaviour becomes stronger, and the ratios are no more constants. 

To sum up, we have shown that Efimov physics for A > 3 manifests with the 
existence of two states for A = 4, 5, 6, and that the ratio between each pair tends to 
be constant as a function of the Efimov angle £. 
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Fig. 1 (color online). Energies of the A = 3 — 6 ground and excited states, E° A ,E^, as a 
function of a -1 , using the two-body gaussian potential (upper panel), and using the two-body 
plus the hypercentral three-body force (lower panel). In both panels we also give the two-body 
ground-state energy E2 calculated with the LM2M2 potential. 
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Fig. 2 (color online). Square root of the ratios between energies of A- particle systems as a 
function of Efimov angle £. The shaded region corresponds to positive scattering length. 
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